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Within the unified model of Bohr and Mottelson we derive the following linear energy weighted 
■ sum rule for low energy orbital 1 + excitations in even-even deformed nuclei 

SS(Mr') S (6/5)e(B(£2;0+ -» 2+ K = 0)/Ze 2 < r 2 > 2 )^ 

with B(E2) the El strength for the transition from the ground state to the first excited state in 
the ground state rotational band, < r 2 > the charge r.m.s. radius squared and e the binding energy 
per nucleon in the nuclear ground state. It is shown that this energy weighted sum rule is in good 
agreement with available experimental data. The sum rule is derived using a simple ansatz for the 
intrinsic ground state wave function that predicts also high energy 1 + strength at 2hw carrying 50% 
of the total mi moment of the orbital Ml operator. 
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Over the last years much experimental and theoretical work has been devoted to the study of orbital 1 + exci- 
tations PJ. Thanks to the combined analyses of (e,e'), (7,7') and (p,p') experiments the occurrence of low lying 
orbital 1^ excitations in even-even deformed nuclei is now a well documented and established fact. Furthermore it 
has been recently found that the summed strength of the observed excitations, that typically appear concentrated 
around an average excitation energy ~ 3 MeV, varies quadratically with quadruple deformation []2|. This last obser- 
vation has motivated theoretical work to understand the origin of the quantitative relation between the orbital 
r* ■ Ml strength and the quadruple deformation parameter j3 from microscopic calculations ||] [|] and to derive general 
formulas connecting Ml and E2 strengths from simplified models 

In Ref. [3] results were presented for summed orbital 1 + strengths in Sin and Nd isotopes that had been obtained 
from deformed HF + BCS calculations with SKIII interaction, using the angular momentum projection formalism 
and substracting spurious contributions ||,[l(J . The results obtained including all possible 2 quasi-particle excitations 
up to E x = 4MeV were found to be proportional to (3 2 , to agree with experimental data, and to account for about 
50% of the total summed strength (obtained including all possible 2 quasi-particle excitations up to E x ~ 25MeV) ||. 
It was also found that for the deformed isotopes the stronger contributions occur in the range 2MeV < E x < AMeV 
and that at higher excitations the strength is much more fragmented. Because calculations neglecting pairing lead to 
summed strengths proportional to (3 rather than to 1 it was argued that dependence on 1 is due to the combined 
effect of pairing and deformation || . Similar results and conclusions have been reached from QRPA calculations based 
on deformed Woods-Saxon potentials with pairing, quadrupole-quadrupole and spin-spin residual interactions . 
The main difference is that the above mentioned QRPA calculations find less orbital 1 + strengths in the high energy 
region. We shall come back to this point latter on. 

In an attempt to establish a general connection between orbital Ml and E2 strengths, linear energy weighted sum 
rules (LEW SR) for orbital 1 + excitations have been recently discussed by Heyde and De Coester jq] and by Zamick 
and Zheng pJ . Within the framework of the interacting boson model (IBM-2) Heyde and De Coester obtain 

X] B(M1- 0+ - lpE x (l+) = ^ B(E2; 0+ - 2+) (1) 
/ / 

with an effective E2 charge appropriate for Boson models. 
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On the other hand Zamick and Zheng, using a quadrupole-quadrupole interaction, find that the LEWSR for orbital 
Ml is proportional to the difference of summed isoscalar and summed isovector E2 strengths. In this work we shall 
show that this approach can be carried further. This will be done by introducing the mean field into the picture. We 
will for the first time look at both the low energy (AN = 0) and the high energy (AN = 2) contributions to LEWSR. 

We shall here consider two closely related models for calculating LEWSR. First a two body quadrupole-quadrupole 
interaction is used to evaluate the double commutator 

S lew (Ml° rb ) = A - Eg . s )\ < f\fl orb \g.s > 

(2) 

= ±tf f <g.s.\[F rb , [H,ft orb ]]\g.s. > 

We note that the commutator in Eq. (2) vanishes for a pure pairing interaction between like nucleons (i.e. between 
protons only and neutrons only) and for a pure spin-spin interaction. Therefore the LEWSR, Eq. (2), should not 
change when said interactions are explicitly considered. The strength of the quadrupole-quadrupole interaction is 
obtained from a self consistency condition. In the second method one evaluates this same double commutator using 
the one body deformed field. 

II. LINEAR ENERGY WEIGHTED SUM RULE WITH A QUADRUPOLE-QUADRUPOLE 

INTERACTION 

We write the quadrupole-quadrupole interaction as 

h q = -xV55> 2 y 2 Mr 2 r 2 )/ (3) 

i>j 

As shown by Zheng and Zamick 0, the value of LEWSR with this interaction is (see also Ref. [16]) 

Si™ (Ml) 



(4) 

= dbxE B ( E2 '> M) - £ B(E2; 1, -1)] 
where B(E2; e„,e n ) is the value of B(E2) for the hypothetical operator 



r 2 Y 2 + e„ r2 Y 2 



The transition is from the J = + ground state to excited J — 2 + states. Eq. (4) holds if we add to the Q.Q 
interaction a pairing interaction between like nucleons. 

For classification purposes it should be noted that for an N — Z nucleus with a T = ground state the value of 
B(E2, 1,1) is four times that of the true electric operator B(E2, 1, 0) for T = final states and is zero for T = 1 final 
states. On the other hand B(E2, 1, —1) will reach only T = 1 final states and is four times the value of B(E2) for the 
electric probe with e p = 1 and e„ = 0. 

A renormalized value of %, which we call xr can be determined, as discussed by Bes and Sorensen lljj ]. The 
interaction of a single valence particle with all the other particles, assuming they have an axially symmetric density 
distribution with the Z axis as a symmetry axis is 

5 

Hdf = -~XRT^-QM(core)(2z 2 - x 2 - y 2 ) (5) 

where Qm(cotc) is the mass Quadrupole moment of the core. It is further shown that the quadrupole moment of the 
core is equal to that of the valence nucleons. Thus the intrinsic mass quadrupole moment Qq 1 is equal to 

Qo 1 = Q M (core) + Q(valence) 2Q M (core) (6) 

If in contrast to the above we use a deformed oscillator hamiltonian 
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(7) 
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then we can make it look like the above expression by introducing the oscillator deformation parameter So 
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(8) 
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where the deformed field term is 



2/ 2 )/3 



We define a deformation parameter S 



5=- 

4 < B? >„ 



(9) 



(10) 



that can be experimentally determined from the measured charge quadrupole moment and r.m.s. radius with bare 
nucleon charges, and assume equal deformation for neutrons (v) and protons (jr). By equating the two expressions 
for Hdf we obtain the following expression for \R 
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Xr 
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muj 2 4ir 
5<R 2 > n 



So/S 



(11) 



However, if we use the QQ interaction amongst all the nucleons in the nucleus, i.e. if we allow AN=2 mixing, then 
the value of x to be used is Xr/^i as a l so discussed by Bes and Sorensen [Q (see pages 143-144). 

We further note that B(E2,1,1) to an isoscalar state is four times that for an electric probe [e p = 1, e n = 0). 



Likewise with e z 
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-1 we reach the isovector states with four times the rate of what an electric probe would 



give. Thus, in order to make easier comparisons with electric probes we give the sum rule in the following way 



S lew (Ml) 



all particles 



9 rnuj S 
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Y J B(E2,l-l) 



(12) 



One final remark to this section. If we had taken for our interaction ^Y^ijQ(i) ■ Q(j), thus keeping i = j terms, we 
would reach the SU(3) limit. The i = j terms add a single particle term of the form r 4 to the potential and cause 
a splitting of single particle energies with different I. In the SU(3) limit in a calculation involving only one major 
shell (i.e. no AiV = 2 mixing) the E2 operator with e p = 1, e n = 1 reaches only one 2 + state-this can be identified 
as the J = 2 + K = member of the ground state band. The E2 operator with e p = 1, e„ = — 1 (isovector) connects 
J = 2 + K = 1 member of the scissors mode rotational band, as was first pointed out by Retamosa et al. |0. It should 
be noted that in shell model calculations the scissors mode band can get fragmented into two bands, one with isospin 
T = To and one with isospin T = To + 1, where To is the g.s. isospin. Thus the isovector E2 will be fragmented to two 
states. Useful formulae for matrix elements of Q.Q in the Boson SU(3) scheme have been obtained by Rosenteel |l^] 



III. THE DOUBLE COMMUTATOR METHOD APPLIED DIRECTLY TO THE DEFORMED FIELD 

We use a simple approximation within the framework of the unified model of Bohr and Mottelson |r| for even-even 
axially symmetric rotational nuclei, where the properties of the ground state rotational band are given in terms of 
properties of the intrinsic ground state. For the intrinsic ground state we use the anisotropic Harmonic Oscillator 
satisfying the selfconsistency condition. The H.O. is used as an auxiliary model because ground state expectation 
values have simple analytical expressions that can then be reinterpreted in terms of those of the "true" intrinsic 
ground state, i.e., in terms of experimentally known properties of the ground state band. 
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Hence in this section we consider a rotational model picture with intrinsic hamiltonian of the form Hdo i n Eq. (7), 
subject to the selfconsistency condition Jl3| 



(13) 



In the rotational model picture the magnetic moment operator is £l = gnl+jlint, with / the total angular momentum 
wihich does not contribute to magnetic excitations and gn the collective gyromagnetic ratio that is experimentally 
determined from the magnetic moment of 2 + , 4 + ... states in the ground state band. In even-even nuclei flint does not 
contribute to static magnetic moments because it has zero spectation value in the intrinsic time even ground state. 
The orbital part of flint is 



rforb 
H-int 



{9i ~ 9r)U = (1 - 9n)L 1T - g R L v 



(14) 



Thus the evaluation of the double commutator in eq. (2) ammounts to the evaluation of the double commutator of 

with Hdo in eqs. (14) and (7), respectively. 
Can we obtain a similar result to that obtained with a Q Q interaction by applying the double commutator method 
directly to the deformed one body field? 
Using the fact that 



we readily find that 



k, \l J ,(2z 2 - x 2 -y 2 ) k 
6(2z 2 - x 2 - y 2 ) k 5i k 5 3 k 

S lew (Ml orb ) = 



(15) 
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(16) 



Using the definition of S in eq. (10) and the relationship for the B(E2) for the ground state K = band of an 
even-even nucleus 



We write eq. (16) in the form 



fl(252;0+->2+* = ()) = —<& 
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S lew (Ml orb ) 



where 
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(17) 



(18) 



(19) 



We compare this result to that of the previous section for the QQ interaction. If we remember that B(E2) W « 
B(E2,l,l)/4), the differences between eqs. (18) and (12) are the presence of the isovector E2 term in the latter and 
also the fact that eq. (18) involves only the B(E2) to the lowest 2 + state and contains the factor F„ u . 

We note that for a nucleus with equal proton and neutron bodies (gn = 1/2, Qq/Qq = 1) the intrinsic orbital 
operator in eq. (14) is purely isovector and the factor defined in eq. (19) takes the value l(F 7ry = 1), as is for instance 
the case for 20 Ne. In this case the dominant contribution to eq. (12) is by far the B{E2) value for the transition to 
the lowest 2 + state, which in the SU(3) limit is the J = 2 + member of the ground state (k = 0) band and exhaust 
the sum rule for quadrupole isoscalar transitions from ground state. 

Thus, in this particular case, one can see that the LEWSR in eq. (18) is similar to that in eq. (12), but in the 
general case one can only check by numerical comparison. One difference is of course that eq. (18) applies to deformed 
nuclei with a well developed rotational g.s. band, while eq. (12) is in principle more general. We would like also to 
point out that the effective nucleon gyromagnetic ratios g\^ — {g\ — gn) for dipole magnetic excitations appearing 
in eq. (14) play a role analogous to the effective charges for dipole electric excitations. In the last case the effective 
charges result from the substraction of the spourious center of mass motion, in the first case g\ ^ results from the 
substruction of the spourious rotation. In what follows we shall see that within the selfconsistent mean field picture, 
implementation of the selfconsistency condition in eq. (13) allows to write eq. (18) in a more practical way for 
phcnomenological analysis. 
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IV. ADDITIONAL CONSIDERATIONS 



First we note that the results in Eqs. (16), (18) can also be obtained considering p.h. contributions. But then we 
get additional insight. If we consider p.h. contributions to the LEWSR 



S lew {Ml orb ) = 

(20) 

( r \ I ^ ^Ulnorb |n \ |2 ,,2 



i;E ph (t P -Ch)\<ph\ f lZ b t.\0>\ 2 



using the selection rules for fl°^ b in the H.O. basis Pq| , we get two sets of p.h. excitations. One set with excitation 
energy \e p — 6h\ = \Tiw x — fiw z \ and one set with excitation energy \e p — = (hw x + hw z ). 
This allows to write Eq. (20) as 

S lew (Ml° rb ) = S LE + S HE (21) 
with Sle the sum from all the low energy p.h. excitations (\e p — eu\ = \hw x — Tiw z \ ~ TiwqSo) 

Sle = [(1 - 9R) 2 Sl E + g 2 R S LE ] /4 (22) 

3 / £ P + I] p \ 

s p le = - nw.) (s^ - - JL ^ JL ) (P + ) 2 (23) 

with p = it,is, and She the sum from all the high energy p.h. excitations (|e p — e/,| = (hw x + hw z ) ~ 2hwo) 

She = [(1 - Sfl) 2 ^ + <?^ £ ] m!v (24) 

^ = ^(^x + hw.) (EJ + ^-^J (P-) 2 (25) 

with 

V V W * 

Using the selfconsistency condition Eq. (13) with £ = it is a simple matter to show that 

87T Jlul 

3 

= ^ ( W^o<2o; p = 7T,v. (26) 



Where we have used the expresion for Qq 



i0 

m \ lo z u>. 
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and eq. (8) for u) x ,uj z . 

Therefore using eqs. (20) to (26) we find the interesting result that 



\s lew {Ml orb ) (28) 
where S lew (Ml orb ) coincides with the value given in Eq. (16), see also eq. (18). 



Sle — She , ) - 
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Thus within this approach the strength function has only two peaks one at low energy (E LE ~ TiWq8q) and one at 
high energy (E HE ~ 2hwo). A pairing interaction between like nucleons, and/or a spin-orbit interaction will cause 
splitting in each of these peaks but will hardly remove strength from the low energy peak to the high energy peak 
and viceversa. Therefore we write a LEWSR for the low energy part of the strength function taking half the value 
given in Eq. (16), or equivalcntly in eq. (18). 

To write the final expression for the energy weighted sum rule for low orbital 1 + excitations we substitute eq. (18) 
into eq. (28). But first we replace the H.O. parameters ujq, Sq as well as Qq/Qq in eqs. (18), (19) by phenomcnological 
parameters directly known from experimental electromagnetic properties of the ground state band. To this end we 
use additional relations that hold within the deformed H.O. model satisfying the selfconsistency condition in eq. (13). 

Contrary to QJ, which is directly related to the experimental quadrupole moments and B(E2) values of the g.s. 
band, there is no direct information on Qq. Using the definition of the collective gyromagnetic ratio gn = I^ r /T cr , 
with X cr = T£ r + I" r the cranking moment of inertia, which in this simplified model is given by 

Or = W4^> (29) 

we can write 



g R = T"jT cr = EVE 



(30) 



and using eqs. (27) and (30) we find that 



QZ 



(1 - 9R)/9R 



On the other hand the charge r.m.s. radius squared in this model is 



< >- 



<Rr>* _ 
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ft fs. 



TiY.' 
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which together with eqs. (8), (10) and (30) give 

5 = 8/ 

and 



mujn 



i S 



9r 



Z <r 2 > D{5) AD{5) Z < r 2 > 



with 

D(5)= (l-ls 
and S fixed by the relations in eqs. (10) and (17), i.e., 

~B(E2;Q+ - 
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2+fc 



0) 



e 2 Z 2 < r 2 > z 



(31) 



(32) 



(33) 



(34) 



(35) 



We note that the experimental quadrupole deformation parameter (3 is usually defined in terms of B(E2) and is 



related to 5 in Eq. (35) by /3 = §£. 

The selfconsistency condition Eq.(13), which ensures that the shape of the potential follows the shape of the density, 
also ensures that i) the momentum distribution is isotropic, and ii) the energy is minimum. These intrinsic ground 
state properties are satisfied by the ground state solution of deformed HF (or HF+BCS) calculations with density 
dependent effective interactions. Thus E^, Wi, can be considered as effective quantum numbers and H.O. frequencies in 
the different directions corresponding to the expectation values of Q Q and r 2 in the "true" intrinsic ground state [|15[ . 



G 



Similarly (|S,Soj/A) can be considered as the effective H.O. energy per particle that corresponds to the true binding 
energy per particle (e — \Eg.s.\/ A), 

hHu/A = 8e/3 (36) 

Substituting Eqs. (31), (34), (36) into Eq. (18), and using Eq. (28), we can finally write a LEWSR for the low 
energy orbital 1 + excitations as 

S le W(M1 or b) _ G 6 Ze B(£2;0+^2+fc = 0) 2 

S LE (M1 )-^u 5D{s) {eZ<r 2 >) 2 »N (37) 

with 

G-ku = g R (l - g R )2A/Z (38) 

For a nucleus with equal proton an neutron bodies G vv = 1 because gn = 1/2 and A/Z = 2. In practice most 
deformed nuclei have gu < 1/2, A/Z > 2 but G^v is still close to 1. For deformed rare earth nuclei since Qr values 
are not known with high accuracy one can use the approximation G^ v ~ 1 to evaluate the right hand side of eq. (37). 
This approximation is used in the results shown in the last columm of Table I. For each nucleus in the table typically 
the value of G^ v oscillates between 0.9 and 1.2, using gn values compatible with the experimental data in Refs. [17] 
and [19]. 

In table I we show results for several deformed nuclei obtained from Eq. (37) using B(E2) values from Ref. [17], 
< r 2 > values from Ref. [18] and binding energies from Ref. [19] approximating G„-„ by one. The results are compared 
to available experimental data on several deformed nuclei. It is interesting to see that this simple approximation leads 
to results in good agreement with experimental data. Therefore Eq. (37) can be considered as a semiempirical 
LEWSR. 

To the extent that more sophisticated microscopic calculations, as those discussed in Refs. [3-5], are able to 
reproduce the low energy part of the strength function for the orbital Ml operator, they also serve as a theoretical 
sound basis to support Eq. (37). As mentioned in the introduction, the main difference between the QRPA results of 
Refs. [4,5] and the results in Ref. [3], that are obtained without inclusion of residual interaction, is that the former 
find less orbital 1 + strength in the high energy region, while the latter find ~ 50% of the scissors mode strength in the 
energy interval AMeV < E x < 2bMeV strongly fragmented in many 2 quasi-particle excitations. On the contrary 
Zawischa and Speth performing QRPA calculations with Migdal interaction find that the strongest scissors like 
excitations take place at E ~ 22MeV . Whether the "true" residual interaction tends to damp down the orbital 1 + 
excitations in the high energy region, or to collect it back into one or a few strong peaks, is still an open question 
that cannot be answered without experimental verification. In this context it should be mentioned that the isovector 
term in Eq. (12), which came from a Q ■ Q interaction, will provide some damping for AN = 2 excitations 0. 



V. CONCLUSIONS 



The main point of this work is to show that with a minimum of number of assumptions one can obtain a simple 
expression (Eq. (37)) which relates the energy weighted sum rule for orbital magnetic dipole strength to the electric 
quadrupole strength, and hence to the deformation parameter 5. The simple expression involves well known quantities 
- the binding energy per particle, the mean square radius and the B(E2) and g^. Although much work along these lines 
has been done using the interacting boson approximation LB. A. we find that we can get equally simple expressions 
working directly with fermions . This indicates that the connection between the orbital Ml and B(E2) is quite general 
and quite natural. 

We further emphasize that the energy weighted orbital magnetic strength should have a low energy part, to which 
our formula, as seen in Table 1, gives very good fits and a high energy part (AN — 2). In this rotational model 
approach the relative energy weighted strengths are easy to calculate. 

Our derivation using the one body field approach may not be as rigorous as the derivation using the two-body 
interaction and it is only applicable to deformed nuclei with a rotational ground state band. The two body interaction 
approach involving an interaction \Q ' Q with \ chosen selfconsistently gives a slightly different expression even for 
N = Z nuclei. In the two body approach we get the difference between summed isoscalar and isovector B(E2)'s. It 
should however be noted that with a Q ■ Q interaction the isoscalar strength goes to only one state in the A7V = SU 
(3) limit - the J=2 + member of the K=0 ground state band. Thus if we drop the isovector term the two approaches 
give the same answer. The isovector term in the AN = SU (3) limit consists of an E2 transition to the K=l, 
J=2 + member of the scissors mode band. It would be of interest to try to measure this B(E2) via electro excitation. 
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However, calculations [[7| indicate that the isovector B(E2) is much weaker than the isoscalar one and so indeed the 
main formula of this work (Eq. (37)) is quite good. 

In the future it will be of interest to extend these energy weighted sum rule techniques to states of higher multi- 
polarity. It would also be of interest to extend the double commutator method from schematic to realistic two body 
interactions. 
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TABLE I. Comparison of experimental values for Y] x E x B(0 + — > 1+) with the results obtained 
from Eq. (37) for the low energy linear energy weighted sum rule of orbital 1 + strength (see text). 
Also shown in the second column are the experimental B(E2) values from Ref. [17] used in this 
work. 



Nucleus 


B(E2) T 


Y,E x B(0+ -> l+)(/&MeV) 


S l £^{Ml orb )(fx 2 N MeV) 




(be) 2 


Expt: E x < 4 MeV 


Theory 


i4b Nd 
148 Nd 
150 Nd 


0.76 
1.38 
2.75 


2.005 a 
3.331 a 
5.95 a ' b 


1.968 
3.339 
6.439 


148 Sm 
150 Sm 
152 Sm 
154 Sm 


0.72 
1.35 
3.44 
4.36 


1.566 c 
3.085 c 
7.003 c 
8.189 c 


1.719 
3.086 
7.310 

8.968 


156 Gd 
158 Gd 
160 Gd 


4.64 
5.02 
5.25 


8.039 d 
8.287 d 
7.246 d 


8.946 
9.616 
9.819 


162 Dy 
164 Dy 


5.28 
5.60 


8.942 b 
9.678 c 


9.926 
10.462 


a Ref. [21]. 
b Ref. [22]. 
c Ref. [2]. 
d Ref. [3]. 



e Ref. [24]. The value quoted has been obtained substructing the spin strength below 4 MeV seen 
in (p,p') and reported in ref. [25]. 
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